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Abstract 

We study kinematic effects due to the approximation of on-shell, collinear partons in shower Monte Carlo event 
generators. We observe that the collinearity approximation, combined with the requirements of energy-momentum 
conservation, gives rise to a kinematic shift, event by event, in longitudinal momentum distributions. We present 
numerical results in the case of jet and heavy flavor production processes measured at the LHC. 
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1. Introduction 

Phenomenological analyses of complex final states pro- 
duced by hard processes at the Large Hadron Collider 
(LHC) rely on event simulation by parton shower Monte 
Carlo generators QJJ] . These are used both to sup- 
plement finite-order perturbative calculations with all- 
order, leading-logarithmic QCD radiative terms and to 
incorporate nonperturbative effects from hadronization, 
multiple parton interactions, underlying events Si. 

The shower Monte Carlo generators S i treat QCD 
multi -parton radiation within collinear ordering approx- 
imations. These approximations have proved to be very 
successful for Monte Carlo simulation of final states at 
LEP and Tevatron. On the other hand at the LHC, unlike 
previous collider experiments, the phase space opens up 
for events with multiple hard scales and multiple jets to 
occur with sizeable rates, while the angular and rapidity 
coverage of detectors extends over a much wider range. 
In this case corrections to collinear ordering can signif- 
icantly affect the structure of multi-jet final states 
i. This, for instance, will influence uncertainties 
of NLO-matched shower calculations S] for jet observ- 
ables. 

In this paper we investigate effects of kinematic ori- 
gin arising from collinearity approximations in the par- 
ton showering algorithms. While the dynamical correc- 
tions studied in |4j-|6|] come from terms beyond NLO 
in QCD perturbation theory (possibly enhanced in cer- 
tain regions of phase space), the contributions studied 
in this paper are not obviously suppressed by powers of 
the strong coupling. Rather, they correspond to approx- 
imations made by showering algorithms in the parton 
kinematics. They can be discussed already at the level 



of leading-order and next-to-leading-order Ollg] shower 
calculations. 

We find that the collinearity approximation, com- 
bined with the requirements of energy-momentum con- 
servation, gives rise to a kinematic shift, event by event, 
in longitudinal momentum distributions. The size of 
this shift depends on the observable and on the phase 
space region, but becomes in general non-negligible 
with increasing rapidities. In Sec. 2 we discuss the phys- 
ical origin of this effect. In Sec. 3 we present numerical 
illustrations for jet production and ^-flavor production 
based on NLO event generators. We give concluding 
remarks in Sec. 4. 



2. Longitudinal momentum shifts in parton showers 

Consider inclusive jet hadro-production. LHC experi- 
ments have measured one-jet cross sections over 
a kinematic range in transverse momentum and rapidity 
much larger than in any previous collider experiment. 
Two kinds of comparisons of experimental data with 
standard model theoretical predictions have been car- 
ried out. The first is based on NLO calculations Hill 
supplemented by nonperturbative (NP) corrections es- 
timated from shower Monte Carlo generators S [loll . 
This shows that the NLO calculation agrees with data at 
central rapidities, while increasing deviations are seen 
with increasing rapidity at large transverse momentum 

A second comparison [9] is based on Powheg calcula- 
tions lfl2ll . in which NLO matrix elements are matched 
with parton showers 

Si 

. This data comparison shows 
large differences in the high rapidity region between 
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results obtained by interfacing Powheg with different 
shower models J2,|3t] and different model tunes lfl3ll . 

The dependence of high-rapidity jet distributions on 
parton showering effects is discussed in 11 4] based on 
results [6] from high-energy factorization jl5ll . These 
results are valid to single-logarithmic accuracy in the 
jet rapidity and the jet transverse momentum, and re- 
sum terms to all orders in the strong coupling a s . In the 
approach J^, [l4ll forward jet production is dominated 
by the scattering of a highly off-shell, low-x parton off 
a nearly on-shell, high-x parton. As noted in IHl . this 
leads to a sizeable fraction of reconstructed jets receiv- 
ing contribution from final state partons produced by 
showering rather than just partons from hard matrix el- 
ements. 

Let us now consider the NLO-matched shower Monte 
Carlo calculations in the light of this physical picture. 
The Monte Carlo first generates hard subprocess events 
with full four-momentum assignments for the external 
lines. In particular, the momenta kj (j = 1,2) of the 
partons initiating the hard scatter are on shell, and are 
taken to be fully collinear with the incoming state mo- 
menta pj, 



kf = XjPj 



0'= 1-2) 
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Next the showering algorithm is applied, and com- 
plete final states are generated including additional 
QCD radiation from the initial-state and final-state par- 
ton cascades. As a result of QCD showering, the mo- 
menta kj are no longer exactly collinear, 



k i * x iPj U = !> 2 ) 
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This corresponds to the soft parton being far off shell 

J~~L 

in the approach On the other hand, in the NLO 
shower calculation the transverse momentum carried by 
partons with kj is to be compensated by a change in 
the kinematics of the hard scattering subprocess. By 
energy-momentum conservation, this implies a reshuf- 
fling, event by event, in the longitudinal momentum 
fractions x, of the partons scattering off each other in 
the hard subprocess. The size of the shift in Xj depends 
on the emitted transverse momenta. 

In the next section we illustrate the longitudinal mo- 
mentum shifts numerically for jet and heavy flavor pro- 
duction using Powheg. 

3. Numerical results 

We consider jet production in the kinematic range J^l] 
and focus on the high rapidity region. In the top panels 



of Fig.Q]we plot the distributions in Xj from Powheg be- 
fore parton showering and after parton showering for the 
low-x and high-x partons. In the bottom panels of Fig.Q] 
we plot the corresponding distributions in the transverse 
momentum k,. 
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Figure 1 : Distributions in the parton longitudinal momentum fraction x (up- 
per plots) and transverse momentum k, (lower plots), before and after parton 
showering, for inclusive jet production at high rapidity. 



We see from the upper plots in Fig. Q]that the kine- 
matical reshuffling in the longitudinal momentum frac- 
tion is negligible for high-x partons but it is significant 
for low-x partons. This effect characterizes the highly 
asymmetric parton kinematics [fl4h . In this region the 
kinematical shift affects the perturbative weight for each 
event as well as the evaluation of the parton distribution 
functions. It will be of interest to examine the impact of 
this phase space region on total cross sections as well. 

The lower plots in Fig.[T]illustrate the distribution of 
transverse momenta in the initial state as a result of par- 
ton showering. We see that this distribution falls off 
rapidly for the low-k, initial state parton, while the dis- 
tribution for the high-A:, initial state parton has a signif- 
icant large-fc, tail. This provides an illustration, at the 
next-to-leading order, of the physical picture described 
in Sec. 2 based on the single-logarithmic, resummed re- 
sults of fi H • 

Analogous kinematic effects can be examined in the 
case of bottom-flavor jet production. The LHC mea- 
surements of Z?-jets llq 1 1 '711 are reasonably described 
by NLO-matched shower generators Mc@nlo 111 80 and 
Powheg lfl9ll at central rapidities, and they are below 
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these predictions at large rapidity and large pj. In 
Fig. [2] we consider /b-jets for several different rapidity 
intervals I1 1611 and plot the gluon x distribution from 
Powheg before parton showering and after parton show- 
ering. We observe similar shift in longitudinal momen- 
tum with increasing rapidity as in the inclusive jet case. 
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Figure 2: Production of /?-jets: distribution in the parton longitudinal momen- 
tum fraction x, before and after parton showering, for different rapidity regions. 



The longitudinal momentum shifts computed in this 
section may affect parton distribution functions in par- 
ton shower calculations. Although we have illustrated 
the shifts using the event generator Powheg, the effect 
is common to other calculation methods using collinear 
shower algorithms. Further studies of this effect will be 
reported elsewhere. 

4. Conclusions 

In this paper we have observed that, in particular phase 
space regions for which inclusive jets and /b-jets are 
measured at the LHC, QCD parton showers are sensi- 
tive to kinematic corrections associated with approxi- 
mations of on-shellness and collinearity on the partonic 
states to which branching algorithms are applied. 

Kinematic reshuffling in the longitudinal momentum 
fractions comes from the implementation of energy- 
momentum conservation in collinearly-ordered shower- 
ing algorithms. It depends on the emitted transverse mo- 



menta, and is present at leading order as well as at next- 
to-leading order. We have illustrated this effect quanti- 
tatively using NLO-matched event generators. 

While the effect is generally small for central pro- 
duction if the center-of-mass energy is not too high, we 
have found that it becomes sizeable in kinematic regions 
probed by hard production processes at the LHC, for in- 
stance jet and heavy-flavor production at high rapidity. 
These regions have not been accessed experimentally 
before at previous colliders (e.g., at the Tevatron). 

The event-by-event shift in longitudinal momentum 
distributions affects both perturbative weights and par- 
ton distribution functions. It contributes in particular to 
the uncertainty on predictions which incorporate results 
of collinearly-ordered showering algorithms. 

Formulations that keep track of non-collinear (i.e., 
transverse and/or anti-collinear) momentum compo- 
nents using unintegrated parton distributions Il20l422ll 
will be helpful to take account of the kinematic effect of 
longitudinal momentum shifts. It will also be of interest 
to investigate the relative size of the kinematic contri- 
bution studied in this article with respect to high-energy 
dynamical fl-Sl effects on jet final states which can be 
included by unintegrated formalisms. 
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